Multifractality of complex networks 
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We demonstrate analytically and numerically the possibility that the fractal property of a scale- 
free network cannot be characterized by a unique fractal dimension and the network takes a mul- 
tifractal structure. It is found that the mass exponents T(g) for several deterministic, stochastic, 
and real-world fractal scale-free networks are nonlinear functions of q, which implies that structural 
measures of these networks obey the multifractal scaling. In addition, we give a general expression 
of T{q) for some class of fractal scale- free networks by a mean- field approximation. The multifrac- 
tal property of network structures is a consequence of large fluctuations of local node density in 
scale-free networks. 



I. INTRODUCTION 

Inspired by the pioneering work by Song et al. 
the fractal property of cornplex networks has been ex- 
tensively studied recently [3-[lOl|. The fractal property 
of a network is measured by the box-covering method in 
which the minimum number of subgraphs (boxes) of di- 
ameter I (in the sense of network distance) required to 
cover the fractal network is proportional to /^■^'' with the 
fractal dimension _Df . Most of real- world fractal networks 
are inhomogeneous in the sense of the scale-free property 
defined by a power-law degree distribution P{k) oc k~'^ , 
where k is the number of connections of a node (degree) 
prj . Thus, the number of nodes in a subgraph of size I 
depends strongly on whether the subgraph includes hubs 
and their neighboring nodes or not, which implies that 
the distribution of local node density is highly inhomo- 
geneous. An inhomogeneous distribution of a physical 
quantity on a fractal object often exhibits the multifrac- 
tal property [T2| - [l5j . In many of fractal objects embed- 
ded in Euclidean space, however, the underlying struc- 
ture seldom shows the multifractal nature because of an 
exponentially thin tail of the mass distribution. On the 
contrary, we expect the multifractal scaling in a scale- free 
network due to large fluctuations of local node density. 
In this paper, we show analytically and numerically that 
fractal scale-free networks (FSFNs) can have the multi- 
fractal property in their structural features. 



II. MULTIFRACTAL ANALYSIS OF 
NETWORKS 

In order to explicate the possibility that a FSFN takes 
a multifractal structure, let us consider, at first, why con- 
ventional fractal structures do not possess the multifrac- 
tal property. In the multifractal analysis, the behavior 
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of a coarse-grained physical quantity on a fractal object 
is argued. Let Xi and /Xi be a physical quantity on the 
discretized position i and its normalized value (measure) , 
this is, 
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The coarse-grained measure (box measure) is then 
given by = X]ie6(/) A^ii where b{l) is a box of size I 
in the system. If the spatial distribution of the physical 
quantity x does not bring a characteristic length scale, 
the qih. moment of the box measure 



(2) 



has a power-law I dependence, namely, (/i^) oc 
where the summation in Eq. ([2]) is taken over boxes of 
size I required to cover minimally the entire system. In 
the case that the exponent T{q) (caUed as the mass expo- 
nent) is a nonlinear (linear) function of q, the distribution 
of the measure is regarded to be multifractal (unifractal) . 
Here, we consider a constant mass of the site i as a physi- 
cal quantity Xi . The normalized measure representing 
the mass density becomes constant and the box measure 
averaged over boxes is proportional to Z^', where 
Df is the fractal dimension of the system. If the fluctu- 
ation of /ifc(i) over boxes is sufficiently small, Eq. ^ is 
approximated as 
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Therefore we have the linear relation r(g) = Df{q — 1) 
representing the unifractal nature of the mass density 
distribution. In fact, a narrow distribution of the box 
measure is widely observed in many fractals embedded in 
Euclidean space p^| - [20j . Thus, most of fractal systems 
take unifractal structures, with some exceptions such as 
mathematical multifractal sets as the two-scale Cantor 
set [2l[ and geochemical distribution of minerals [22| - [2i | . 
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FIG. 1. (Color online) Two types of covering schemes for 
the (w, v) flower with u — v = 2 in the third generation, 
(a) The network is covered by 16 {u, v) flowers in the first 
generation (scheme I), and (b) the network is first covered 
by four subgraphs of size / — 2 centered at the largest hubs 
(white nodes) then covered by eight subgraphs centered at 
the second largest hubs (gray nodes) (scheme II). 



The box measure of node density in a FSFN, how- 
ever, has large fluctuations over boxes [l3|. This is be- 
cause fJ.b{i) of a box including a hub node can be much 
larger than that of a box without hubs. If the proba- 
bility distribution function of the box measure has a fat 
tail like a power-law or log-normal form, cannot be 
approximated by {^J,b(l)) and the failure of Eq. ^ gives 
a possibility of the multifractal scaling in structures of 
complex networks. 

When applying the multifractal analysis to complex 
networks, it is unavoidable that covering boxes (sub- 
graphs) of network diameter I overlap each other |5|, if 
every box includes the maximum number of nodes un- 
less the subgraph diameter does not exceed I as in the 
conventional multifractal analysis [l^. Therefore, the 
measure fii defined by Eq. ([1]) cannot be normalized as 
J2b{i) t^b(i) ~ ^^^"^ mass exponent does not satisfy 
the basic condition r(l) = 0. In order to overcome this 
difficulty, we modify the definition of the measure as 



In this definition, the normalization constant, and then 
/ii, varies with the box size I. It is easy to confirm that 
the mass exponent calculated from Eq. (|4]) satisfies the 
general conditions t(0) = —D[ and t(1) = and for a 
unifractal system T{q) = Df{q — 1). We set hereafter 
Xi = 1 to analyze the multifractal nature of the node 
density in a complex network. 

III. MULTIFRACTALITY OF {u, v) FLOWERS 

Using the above definition of fii, we first examine an- 
alytically the multifractal property of the (w, w)-fiower 
model [9|. This deterministic model provides a class of 



FSFNs. In this model, we start with the cycle graph con- 
sisting of w = u + V (1 < u < v) nodes and edges (the 
first generation). The network in the nth generation is 
obtained by replacing each edge in the (n — l)th genera- 
tion network by two parallel paths of u and v edges. The 
network with large n has the scale-freeproperty with the 
degree exponent 7 = 1-1- log w/ log 2 [9|. The number of 
nodes in the nth generation network is given by 

and the number of edges is w". In addition, when w is 
even, the diameter is written as 

Ln = + ' . 6 

\2 u — 1 J u — 1 

From the network formation algorithm, the nth genera- 
tion network is constructed by w""™ of (u, v) flowers in 
the mth generation, where m < n. This implies that if we 
cover the nth generation network by the mth generation 
subgraphs as shown in Fig. 1(a) the number of covering 
boxes iVh(L„) is Using Eq. dH), we can rewrite this 

relation as A^b(L„) = w"[(L„ -f 6)/a]-^°s'"/i°s", where 
a = {w/2 + b)/u and b = {v~u)/{u~l). If the length Lm 
is large enough, we have the fractal dimension given by 
D{ = log w/ logu The above covering scheme (named 
as the covering scheme I), however, does not lead the min- 
imum value of Ni,(L^^ for a fixed Lm- Let us consider the 
following covering scheme (named as the covering scheme 
II). At first, we cover the {u,v) flower by subgraphs of 
size Lm centered at the largest hubs, then the remaining 
network is covered by subgraphs centered at the second 
largest hubs. Repeating this procedure until all nodes are 
covered by subgraphs [as shown in Fig. 1(b)] the number 
of covering boxes Nb[L^) becomes less than w"^™. Al- 
though two covering schemes yield different Nb[L,n)^ these 
are in proportion as shown later, and then the fractal di- 
mension calculated by the covering scheme II is the same 
as that by scheme I. Therefore, both covering schemes are 
valid for the fractal analysis. In the multifractal analysis, 
however, we treat not only Nb(^i^^^{= (/i^"")) but also 
(^^ ) for any q. Since the moment (/x^ ) calculated 
by the covering scheme I is generally not proportional to 
) by scheme II, we need to choose scheme II with less 
(probably minimum) covering boxes for the multifractal 
analysis of the (u, v) flower. It is, in general, quite im- 
portant to cover a network by less number of boxes in the 
multifractal analysis comparing to the case of the fractal 
analysis. 

Let us cover the {u, v) fiower in the nth generation 
by boxes of size I = Lm (1 m ^ n) in the covering 
scheme II. The number of boxes iVf,(s_i^) centered at the 
sth largest hubs is equal to the number of such hubs, thus 
we have 

Nb{s,Lr„) =Vs- Vs-1 {1 < s <n~m) , (7) 
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where b{s,Ljn) represents a box of size centered at 
the sth largest hub and ~ 0. Since the number of 
nodes i>s{Lm) in b{s, Lm) is presented by 

i>,(i„) = 2"-™-^+i(^.„/2 + 1) ^ 2"-"-^™ , (8) 

the total number of nodes in all boxes, namely, the de- 
nominator of Eq. gl), is given by Nh{s,L„,)Vs{Lm)- 
Using Eqs. JS]), ([7]), and ([5]), the measure defined by 
Eq. Q with Xi = 1 is thus expressed as 

ui — 1 
" u;"(i«-2) • 

It should be noted that Hi for the (u, v) flower is inde- 
pendent of the box size though the box measure given 
by Eq. ^ generally depends on I. Denoting //^ inde- 
pendent of i by /i, the qth moment of the box measure 

is calculated from {^ilj = Eh(s,L„o (Eiefc(s,L„) ^^ = 

YTsZT ^b(s,L^)[^^Vs{L,-a)\'i. By means of Eqs. (P and 
([7]) — (O, the quantity ) is calculated as 



(Mi, 



w 



n(l-g) 



L.m. + b\ 



(10) 



where a and h are defined below Eq. ([6]), W^g = (w - 
2'^)/ {w — 2), and the approximation ~ w™{w—2)/{w- 
1) is used. The dominant term in Eq. ()10p for large L^i 
depends on q, and the mass exponent is given by 



log(w/2) . logw 
q — 11 <Z > 



r{q) 



(9-1) 



logM 

logui 



logw 



if g < 



log 2' 
log w 
bi2' 
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The nonlinear form of r(g) indicates that local node den- 
sities of the (m, v) flower are distributed in a multifractal 
manner (25j . Although the above argument holds only 
for even values of w, we found that Eq. pT|) is a good 
approximation also for odd w. It should be noted that 



(A^iJ 



w 



"(1-9) 



[{Lm + 6)/a](9-i)i°s»/iogt^ obtained by 



scheme I is proportional to Eq. (fTO|) for q < \ogw/ log 2, 
but not otherwise. 

It is generally difficult to find the way to cover mini- 
mally a given network because the minimization of the 
number of covering boxes is known to be NP hard. We 
then need to cover the network by a less number of boxes, 
as an approximation, in actual multifractal analyses. A 
variety of such covering methods have been proposed so 
far [l[ [^-Q. In order to clarify whether such covering 
techniques proven to be efficient for fractal analyses still 
work even in multifractal analyses sensitive to the cover- 
ing way, we compare the analytical expression Eq. pT|) 
with the numerically calculated T(q) by adopting the 
compact-box-burning algorithm [4] modified to shorten 
the computing time 0- The results shown in Fig. 2(a) 
clearly demonstrate the validity of this covering method. 
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FIG. 2. (Color online) Mass exponent r(g) for several fractal 
complex networks, (a) Solid line (red line) and circles indi- 
cate r(g) for the {u,v) flower with u — v — 2 calculated by 
Eq. (jlip and by the numerical box-covering for the network 
in the eighth generation (A'^ = 43,692), respectively. Squares 
represent numerical results of r(g) for the network formed 
by the minimal model with m = 2 in the seventh generation 
(TV = 15,626). Dashed line (blue line) indicates the theo- 
retical curve for this minimal model predicted by Eq. (I13|l . 
Triangles and crosses show r(g) for the giant components of 
the fitness model {N = 100, 000) and the Erdos-Renyi ran- 
dom graph (A'^ = 200, 000) at their percolation transitions, 
respectively. The degree exponent for the fitness model is set 
to be 7 = 4.0. We averaged (/i^) over 100 realizations both for 
the fitness and the random graph models. Dashed-dotted line 
(green line) indicates r(g) given by Eq. (|13|l with 7 = 4.0 and 
D{ = 2 corresponding to the fitness model. Results except 
for the {u, v) fiower are shifted vertically for clarity though 
all T(g) actually pass through t(1) = 0. (b) The mass expo- 
nent T{q) for the WWW is presented by circles. Straight lines 
indicating the slopes of r(q) for q <C — 1 and q ^ 1 are just 
guides to the eye. 



IV. MEAN-FIELD ARGUMENT 

Let us generalize our argument to some extent. We 
treat a FSFN of nodes, whose degree distribution is 
given by P{k) cx k~'^ . Covering the network minimally 
by iVb(i) boxes of diameter Z, the mean number of nodes 
{vi) in a box is given by N/Nm-^, where -/Vh(;) oc 
Regarding mutually connected boxes of size I as a renor- 
malized network ,lj], the fractal property of the original 
network assures the relation Pi{ki) oc k'j^'^ , where Pi{ki) 
is the degree distribution function of the renormalized 
network. Here, we assume that each renormalized node 
is statistically equivalent and the number of nodes vi^ki) 
in a covering box corresponding to the renormalized node 
of degree ki has negligibly small fluctuations over boxes. 
In this mean- field approach, the quantity vi^ki) is repre- 
sented by its mean value. 



(12) 
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where (fc) = (ki) is the average degree of the original 
network. The {u, v) flower satisfies Eq. rigorously 
in the thermodynamic limit (n — >■ oo). Since the box 
measure /if,(i) is given by vi{ki) normalized by N, this is, 
Mb(z) = ki/[{k)Nb(i)] oc {ki/{k))l°f, the qth moment (^^) 
can be calculated by Eq. ([2]). Considering that the max- 
imum degree fcmax is proportional to N^/'--''~^\ we have 
cx for g < 7 - 1 and (^«) oc Z'3Cf(7-2)/(7-i) 

for g > 7 — 1. Therefore, the mass exponent r(g) of this 
network is presented by 

(iq-l)Df if9<7-l 

This result implies that FSFNs satisfying Eq. gen- 
erally take multifractal (bifractal) structures. Our result 
Eq. prj) for the (u, w) flower is a special case of Eq. ([T^. 

In order to confirm the validity of the above mean- 
field argument, we numerically calculate the mass ex- 
ponent for the minimal model proposed by 0], employ- 
ing the compact-box-burning algorithm [J|. A network 
formed by this model possesses the scale-free property 
with the degree exponent 7 = 1-1- log(2TO -I- l)/logm 
and takes a fractal structure with the fractal dimension 
Df = log(2m -I- l)/log3. The minimal model also satis- 
fies Eq. as in the case of the {u,v) fiower model 2]. 
The nonlinear behavior of numerically calculated r(g) in- 
dicated by squares in Fig. 2(a) agrees quite well with the 
theoretical prediction Eq. (IT^ . As an example of net- 
works not satisfying Eq. ([T^ due to large fluctuations of 
iyi{ki), we treat a stochastic model of FSFNs, namely the 
fitness model proposed by [l^ . A network formed by this 
algorithm also has the scale-free property and the giant 
component exhibits the fractal nature at the percolation 
transition 0. Triangles in Fig. 2(a) represent T{q) for 
FSFNs formed by the fitness model with 7 = 4.0 and 
Df = 2 0. The exponent T{q) is also a nonlinear func- 
tion of q, which suggests the multifractal structure of the 



network, but cannot be described by Eq. (IT3l) . Finally, 
we calculate the mass exponent t(o) for the World Wide 
Web (WWW) with 325, 729 nodes [I^, which is known to 
be a representative real- world FSFN fl'|. Although the 
nonlinearity of T{q) is weak as shown in Fig. 2(b) and 
the result is not described by Eq. (IT51) . two tangential 
lines in the extreme regimes q —1 and 5 3> 1 have 
definitely different slopes, which shows the multifractal 
structure of the WWW. The multifractal property found 
in these networks is obviously caused by the scale- free na- 
ture of networks. In fact, r(q) for the giant component 
of the Erdos-Renyi random graph [2^ at the percolation 
threshold is, as shown by crosses in Fig. 2(a), a linear 
function of q, where the giant component is fractal but 
not scale free. 

V. SUMMARY 

In conclusion, we demonstrated analytically and nu- 
merically that fractal scale-free networks (FSFNs) can 
take multifractal structures in which the fractal dimen- 
sion is not enough to characterize fractality of systems. 
The multifractal nature is caused by large fluctuations 
in local node density in scale-free networks. Although 
all examples treated in this paper exhibit the multifrac- 
tal nature, further investigations are needed to clarify 
whether any FSFNs take multifractal structures. It is 
also crucial to study how the multifractal property affects 
physical phenomena or dynamics on complex networks. 
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